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ABSTRACT 

Radiative diffusion damps acoustic modes at large comoving wavenumber (k) before decoupling ("Silk 
damping")- In a simple WKB analysis, neglecting moments of the temperature distribution beyond the 
quadrupole (the tight-coupling limit), damping appears in the acoustic mode as a term of order ik 2 T~ l , where 
f is the scattering rate per unit conformal time. Although the Jeans instability is stabilized on scales smaller 
than the adiabatic Jeans length, I show that the medium is linearly unstable to first order in f" 1 to a slow 
diffusive mode. At large comoving wavenumber, the characteristic growth rate becomes independent of spatial 
scale and constant: (ficH^) -1 ~ (1287rG/9Krc)(p m /pft), where a is the scale factor, p m and pi, are the matter 
and baryon energy density, respectively, and kj is the Thomson opacity. This is the characteristic timescale 
for a fluid parcel to radiate away its total thermal energy content at the Eddington limit, analogous to the 
Kelvin-Helmholz (KH) timescale for a radiation pressure-dominated massive star or the Salpeter timescale 
for black hole growth. Although this mode grows at all times prior to decoupling and on scales smaller than 
roughly the horizon, the growth time is long, about 100 times the age of the universe at decoupling. Thus, 
it modifies the density and temperature perturbations on small scales only at the percent level. The physics 
of this mode in the tight-coupling limit is already accounted for in the popular codes CMBFAST and CAMB, 
but is typically neglected in analytic studies of the growth of primordial perturbations. The goal of this work 
is to clarify the physics of this diffusive instability in the epoch before decoupling, and to emphasize that 
the universe is formally unstable on scales below the horizon, even in the limit of very large f. Analogous 
instabilities that might operate at yet earlier epochs are also mentioned. 

Subject headings: cosmic microwave background — cosmology : theory — large-scale structure of universe 
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1. INTRODUCTION 

In standard analytic treatments of perturbation theory in the 
epoch near radiation-matter equality the tight coupling ap- 
proximation is employed, which amounts to neglecting mo- 
ments of the temperature distribution beyond the quadrupole. 
In this limit, the scattering rate per unit conformal time 
f is very large and, neglecting gravity at large comoving 
wavenumber, the dispersion relation for acoustic modes is 
easily derived; the medium supports stable acoustic oscilla- 
tions modified by radiative ("Silk") damping on small scales 
(Silk 1967, 1968; Peebles & Yu 1970; Weinberg 1971; Ma & 
Bertschinger 1995; Hu & Sugiyama 1995ab, 1996, 1997; Do- 
delson 2003). In the WKB approximation the damping term 
enters as a complex correction to the acoustic mode frequen- 
cies that is cx k 2 /f, becoming increasingly important at small 
spatial scales (see, e.g., Dodelson 2003, §8.4). 

Most analytic calculations of radiative damping of acoustic 
modes neglect gravity at large comoving wavenumber, k, be- 
cause the frequency of acoustic oscillations (w « c s k, where 
c s is the sound speed) is large and the medium is Jeans sta- 
ble. Here, I show that on scales smaller than the Jeans length 
the medium is linearly unstable to a slow diffusive mode that 
is of order f" 1 . Thus, on small scales where the medium is 
dynamically Jeans stable and supports stable (but radiatively 
damped) acoustic oscillations, it is unstable to an orthogo- 
nal diffusive mode. This mode was earlier discussed by Ya- 
mamoto et al. (1998). The physics of this slow diffusive mode 
is already included in CMB codes like CMBFAST and CAMB, 
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which include gravity on small scales in the tight-coupling 
limit in their solution to the coupled equations for perturba- 
tions, and in earlier work employing full numerical solutions 
to the coupled perturbation equations. However, the physics 
of this mode has been often neglected in analytic studies of 
the growth of perturbations, and for this reason a more com- 
plete discussion is warranted; it may affect analytic estimates 
for the characteristic baryon-acoustic oscillation (BAO) scale, 
the physics of low-mass dark matter halos, and the transfer 
function. 

In §12 I review the physics of acoustic modes and radia- 
tive damping in the cosmological context, and I show that the 
unstable diffusive mode appears in a simple WKB treatment 
that includes gravity on small scales. The discussion is re- 
lated to that of the "terminal velocity" mode in Yamamoto et 
al. (1998). I emphasize that the growth of this mode is not tied 
to the breakdown in the coupling between gas and radiation; 
indeed, it appears at the same order in r as Silk damping, and 
a temperature difference between the radiation and the gas is 
not required for growth. I note that analogous instabilities 
should operate in the epoch of neutrino decoupling, or when- 
ever the medium is partially supported against self-gravity by 
a relativistic and diffusing particle species. Section|3]presents 
a complimentary analysis in an expanding background clos- 
est in spirit to Silk (1967), and to the discussion of adiabatic 
modes in standard textbooks. In §H] I provide a brief sum- 
mary. 

2. ACOUSTIC MODES & SILK DAMPING 

Following Hu & Sugiyama (1996), the equations for the 
evolution of the moments of the photon temperature distribu- 



tion 9„ (n = 0, 1 , 2) in the tight-coupling limit are 
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where Vj, is the baryon velocity and the last expression ne- 
glects 02- The continuity and Euler equations for the baryons 
and dark matter are 
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respectively, where R = 3p b /4p r , 

f = dr /di] = -n e (TTa \ 
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rj denotes conformal time, kt is the Thomson opacity, and 
the notation of Dodelson (2003) has been adopted. The field 
equations can be written as 
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and 
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where p b , p m , and p r denote the baryon, matter, and radia- 
tion energy density, respectively. The approximate equality 
in equation (HJ assumes that one consider epochs sufficiently 
before recombination that the universe is fully ionized. The 
above set of equations neglects moments of the temperature 
distribution above the quadrupole, the polarization hierarchy, 
and neutrinos. 

In the simplest derivation of Silk damping, one neglects 
both the gravitational terms and the expansion of the universe 
(so that d/a = $ = $ = "J = 0) as well as the matter. Expanding 
the time dependence as e' ut , one finds that 
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Expanding to lowest order in the limit of large f, there are two 
acoustic modes 
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where c 2 = c 2 /[3(l +R)], and the purely damped mode (see 
e.g., Blaes & Socrates 2003) 2 

l+R\ ik 2 R (8-/?) 
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Because f is always negative and because I have written the 
time dependence as e IUJt , all terms in all three modes involv- 
ing f represent damping. Note that equation (fTTT i contains a 
fourth mode that has uj = 0, which I now discuss when gravity 
is included in the system. 

2 Note that champed disappears from the analysis if the gas and radiation 
temperatures are assumed to be equal, as in {J3](e.g., Blaes & Socrates 2003; 
Kaneko & Morita 2006; Appendix B of Thompson 2008). 



2.1. The Simplest Case with Gravity 

The neglect of gravity in the above derivation is justified on 
small scales for acoustic modes whose frequencies are oc c s k 
(eq. Ifl2l ). It turns out that it is also justified for the purely 
damped radiative mode (eq. fT3l ). However, including gravity 
on small scales introduces a qualitatively different behavior to 
the fourth mode, which is purely damped on scales larger than 
the Jeans length and purely unstable on scales smaller than the 
Jeans length. 

Keeping the field equations, but again ignoring terms oc 
d/a, oc <I>, and the equations for the matter, I find 
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where £ 2 = l6irGa 2 p r . As before, there are four orthogonal 
modes admitted by equation (fPfl) The first and second are 
gravity- and radiative-diffusion-modified acoustic modes. To 
first order in f" 1 these are 
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This expression should be compared with equation ( fl2b . Note 
that taking either the k — > oo or the £ 2 — > limits in equation 
( fTSI l recovers equation ( fl2b . Furthermore, the first term in 
equation (TT3T > is identical to what one would expect for the 
Jeans instability. Setting the first term in square brackets to 
zero, I derive the Jeans length: 
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The third mode admitted by equation ( [Pil l is the strongly 
damped diffusion mode, uidamped, of equation ( fT~3b , which is 
unmodified by the inclusion of gravity to second order in f -1 . 
The fourth and final mode admitted by equation ( fl4l ) is purely 
imaginary. It is 

ik 2 ( R 2 \ ( i 2 . 
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to first order in f . This expression follows from equating 
the last two terms in equation (TT~4-b . Note that in the limit that 
the medium is Jeans stable c 2 s k 2 > (1 +R)£ > 2 /3 (large comov- 
ing wavenumber, small spatial scale) equation ( fTTl ) implies 
that the medium is unstable. Conversely, on scales larger than 
Aj, Lu Kli is damped. In the high-A: limit (c 2 k 2 > (1 +R)£/3), 
equation ( fTTl ) can be written simply as 
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which is independent of spatial scale. Precisely when the 
medium is stable to the classical Jeans instability it is unstable 
to a slow diffusive mode. 
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The characteristic frequency is smaller than that for any of 
the other modes: the timescale for growth is [3itG/(ktc)]~ 1 ps 
2 x 10 16 s, much longer than the radiation acoustic timescale, 
^acoustic ~ ( c sk)~ l , on a scale of order the Hubble radius and 
the interaction timescale damped- The subscript "KH" is 
used to emphasize the connection with the Kelvin-Helmholtz 
timescale for gravitational contraction of a radiation pressure 
supported self-gravitating body (see ^2. 2b . The timescale ?kh 
is roughly m e /m p shorter than the timescale for the radia- 
tive instability derived by Gamow (1949) and Field (1971) 
(~ 10 12 yr; their eqs. 20 & 58, respectively), and does not for- 
mally require a temperature difference between the radiation 
field and the gas (^3). 

2.2. The Simplest Case with Dark Matter 

The slow diffusive mode derived above changes when dark 
matter is included primarily because the strength of the poten- 
tial changes. In the WKB limit, the mode structure becomes 
more complicated when the Euler and continuity equations 
for the dark matter are included and the dispersion relation 
is cumbersome. Nevertheless, neglecting a /a, $ and neutri- 
nos as before, the slow diffusive mode can be isolated in the 
dispersion relation. To first order in f _1 , it reads 

k 2 ( R \ r 4 2 + £ 2 „(l + 8/(9fl)) 



3r V l+R 



[c^-^(i+R)/3-a. (I9i 

where £f mr = l6nGa 2 p bmr . Compare equation ([19} with 
equation fl7\ . Again taking the strongly Jeans stable limit 
of small spatial scales and large k (c 2 k 2 ^> £^+£ 2 (l+i?)/3), I 
have that 
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In the limit of small R this expression becomes 
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As in equation (TTHt . this timescale is independent of spatial 
scale. As long as the medium is diffusive (the tight-coupling 
approximation holds) and one consider scales smaller than the 
Jeans length, the medium is unstable. Furthermore, note that 
taking instead the large R limit in equation (l20l recovers the 
basic form of equation (II St . but with the substitution pi, — ► 
pb + p m in the numerator. However, for z > 1000, R < 1 and 
equation (f2Tb is the relevant limit. The characteristic growth 
timescale is then 
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Note the correspondence with the Salpeter timescale for black 
hole growth. The characteristic timescale in equation (l22b can 
be understood as the timescale for a region of differential vol- 
ume to radiate away its total (relativistic) thermal energy con- 
tent (e ~ pbc 2 ) at the Eddington limit (e ~ 4irGp tot c / kt)'- 3 
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This is simply the Kelvin-Helmholz timescale. f KH can equiv- 
alently be thought of as the ratio of the square of the dy- 
namical timescale, (47rGjOt ot ) _1 , to the photon-baryon inter- 
action timescale, (KrPbc) ■ Alternatively, one may think 

3 The "9" in equations 12 It & i22\ originates in the numerical factor in 
equation {3) connecting 0] and ©2 ( see > e -g-. Hu & Sugiyama 1996). 



of the growth timescale as ?kh = fdiff (%n/*r) , where taaf = 
ck 2 j '(3kt Pb) is the diffusion time, % n = (47rGptot) -1 ^ 2 is the 
dynamical time, and t r = (c s k)~ l is the radiation acoustic 
sound-crossing timescale on a scale k~ l . Finally, one can also 
obtain this characteristic timescale by equating the frictional 
force between the radiation and the fluid due to Thomson scat- 
tering with the gravitational force. The timescale for instabil- 
ity is ~ m e /m p shorter than the timescale for the instability 
discussed in Gamow (1949) (see also Field 1971). Note that 
because the mode that becomes unstable has uj = in the adia- 
batic limit, one may identify this unstable mode as the entropy 
mode (e.g., Lithwick & Goldreich 2001). 

Importantly, ?kh in equation (l22l is roughly 100 times the 
age of the universe at decoupling and approximately 500 
times the age of the universe at radiation-matter equality. At 
earlier times the instability also operates, but because ?kh is 
a constant in the large-A; limit, it becomes increasingly long 
compared to the Hubble time at higher redshift. This is most 
clearly seen by writing 



(24) 



where T is the interaction timescale and H is the Hubble 
parameter. Because T ^> H(z) in the tight-coupling limit, 
tjQiH(z)a is much larger than unity, and the instability is slow. 

For a non-relativistic discussion of this unstable mode and 
its potential astrophysical applications, see Kaneko & Morita 
(2005) and Thompson (2008). The same mode has been dis- 
cussed by Yamamoto et al. (1998), where they term it the "ter- 
minal velocity mode" (their Section 3.3), since one can write 
the baryon velocity as Vj ~ k^R/r (their notation) by equat- 
ing the gravitational force with the drag force associated with 
the coupling to the radiation field. As I have shown above, the 
mode is unstable at arbitrarily large (but, not infinite) f, with 
fixed growth timescale, and operates until decoupling. As Ya- 
mamoto et al. (1998) show, the instability smoothly joins the 
pressure-less collapse solution after decoupling. 

Note that physically equivalent modes should operate 
whenever the medium is supported against self-gravity by a 
particle species that diffuses. For example, in the epoch be- 
fore neutrino-matter decoupling at a time ~1 s and at a tem- 
perature larger than ~ 10 10 K a similar mode is expected to 
have growth rate of order 
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where here p m is the mass density in all matter, and T ~ 
G 2 F T 5 jh is the weak interaction rate; T w 0.27j^ eV s" 1 where 
?MeV = T/MeV. The mass density in matter at z ~ 10 10 is of 
order p m ~ 1 g cm" 3 and so I find that (u>/a) ~ —i 10~ 6 s" 1 . 
That is, the growth timescale is of order one million times the 
age of the universe at that epoch. 

3. SLOW GRAVITATIONAL INSTABILITY 
IN AN EXPANDING BACKGROUND 

The preceding analysis approximated the background state 
as fixed. However, the growth timescale derived exceeds the 
timescale for the background state to change by a factor of 
more than 100. This invalidates the WKB calculation and 
necessitates a treatment with an evolving background. This 
is most simply treated on small scales where the gravita- 
tional potential may be approximated as Newtonian. For sim- 
plicity of presentation, I treat the non-relativistic case with 



aT 4 <C pc 2 and I neglect dark matter. In this case, the Eule- 
rian equation governing the evolution in an inertial frame are 
(e.g., Silk 1967) 
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Here, p is the mass density of gas, u is the velocity, and 
p = aT 4 /3 and e = 3p are the radiation pressure and energy 
density, respectively. The gas and radiation are assumed to 
have the same temperature, and it is assumed that p greatly 
exceeds p g = nk^T ', the gas pressure. 4 

In equations d26t - d30l ). I have adopted the notation of Pee- 
bles (1993) (Chapter 5). The subscript r reflects reference 
to an inertial coordinate system. Transforming to a comoving 
coordinate system withx = r/a(f) and u = a'x+v(x,f) = dx+ax, 
where v is the peculiar velocity, and writing p = p(t)(l + 
5(x,t)), with p(t) oc a(t)~ 3 , and taking d/a = -(A/3)-KGp(i), 
equations (I26ii-(l3"0"il become 
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where time and space derivatives are now understood to be 
in the comoving frame and $ = 4>(x,t) + (2/3)nGp(t)a 2 x 2 . 
In analogy with the density perturbation I write e = e(f)(l + 
e(x,f)), and equation ( l33l can be rewritten as 
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So far the equations written are completely general for the 
problem at hand. I now specialize to the case of small depar- 
tures from the background state so that terms of order 5|v|, 
Ss, and e|v| can be neglected. The resulting equations can be 
written as 
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Note that the left hand side of equation d37l i is simply 

de A 88 _ 4 d /As' 
lh~3~dt~3dt 

4 For the general case with gas pressure, see Thompson (2008). 
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where s is the total entropy. Writing e = Ae/e = 4 AT /T, 
equations d36l l and fl37T i become identical to equations (1) & 
(2) of Silk (1967). Expanding these expressions in Fourier 
components with wavenumber k = k/a, taking a(t ) oc t 2//3 , and 
defining the temperature perturbation as q = e/4, results in two 
coupled ordinary differential equations for the time evolution 
of the density and temperature perturbations: 



and 
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^-3^ + ^ = 0, 
where over-dots denote time derivatives, 

. _ ck 2 
3kp 

is the diffusion rate, 

2 4 e dp 
Cr= 9~p = d~p , 

is the adiabatic radiation pressure sound speed, and the sub- 
script k has been added to emphasize that these equations are 
for the evolution of the density and temperature, 6 and q, at a 
specific comoving wavenumber. 

3.1. Adiabatic Modes 

In the adiabatic limit u> = in equation d40l i. d/dt(As/s) = 
0, q = (1/3)5, and the system becomes (e.g., Peebles 
eq. 5.124) 



(43) 



In the long-wavelength limit (k — > 0) one recovers the classic 
equation for the time evolution of the density contrast, mod- 
eled as a purely pressure-less fluid, which has the growing and 
damped solutions 5 + oc f 2 / 3 and 5- oc t~\ respectively. Taking 
k — > 00, equation d43l describes an acoustic wave, damped by 
the expansion of the universe. 

3.2. Non-Adiabatic Evolution 

The full non-adiabatic evolution is more complicated. First, 
as shown by Silk (1967) and others, the acoustic modes are 
damped by radiative diffusion. Second, — the purpose of this 
paper — as shown in the previous sections, it is formally in- 
correct to neglect gravity at large comoving wavenumber, be- 
cause to do so is to effectively set the driving term for the 
diffusive gravitational instability to zero. 

To illustrate the importance of these effects, I solve equa- 
tions (O and (|40j for & k (t) and qffl with 5\(t = 0) = 0, 
§- k (t = 0) = 1, qjfi = 0) = 1/3, c 2 l 2 = 100/f 2 / 3 , and u = t 2 /VTo 
over the range 0.1 <t < 100. These parameters ensure that 
the mode considered is much smaller than the Jeans scale and 
that the diffusion rate is small. Figure Q] shows results for 
the time evolution of the density and temperature perturba- 
tions with and without gravity (heavy and light solid lines, 
respectively). The acoustic mode is damped by both the ex- 
pansion of the universe and Silk damping. At early times, the 
evolution with and without gravity is very similar. At later 
times, the presence of the unstable growing mode becomes 
more pronounced and the evolution of 5~ k and q~ k is qualita- 
tively different. The right panel shows a zoomed-in version of 
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FIG. 1 . — Left Panel: Solution to equations (39) and )40t at large k as a function of time both with gravity (the "2/3r 2 " term; heavy solid line) and without 
gravity (light solid). The top and bottom panels show the evolution of the density and temperature perturbations, &~ k (t) and qi(f). The acoustic wave is rapidly 
damped by both the expansion of the universe and diffusive damping. At early times the two solutions track each other closely. After the mode is damped, the 
evolution is qualitatively different as a result of the slow diffusive mode. Right Panel: A zoom in on the late-time evolution of 8~ k to highlight the difference 
between the calculations with and without gravity. The increase in & k (t) at late times is not a result of the breakdown of the tight-coupling approximation, but 
is instead a result of the instability highlighted in this paper. Importantly, in a more complete calculation the medium becomes optically-thin, the radiation 
decouples, and the solution joins the classic pressure-less free-fall collapse solution. See Yamamoto et al. (1998) and Singh & Ma (2002) for detailed calculations 
at small scales. 



the late-time evolution of 5~ k . At t « 6, as the mode becomes 
fully damped, the fractional difference in 5~ k between the two 
solutions is a few to 10%. Of course, this example does not 
include the physics of decoupling of the mode and gravita- 
tional free-fall into the dark matter potential, which would be 
included in a full calculation, and so the difference between 
the two solutions becomes exaggerated at late times (see Ya- 
mamoto et al. 1998; Singh & Ma 2002). 

4. DISCUSSION 

The original treatments of acoustic modes near the epoch 
of radiation-matter equality focused on the importance of ra- 
diative damping of acoustic modes at large comoving wave 
number. The general equations derived by Silk (1967) (and 
Silk 1968; Peebles & Yu 1970; Weinberg 1971) contain the 
physics described above, but the slow diffusive mode these 
equations admit was not discussed because gravity was ne- 
glected when calculating the damping of acoustic modes on 
small scales. This approximation eliminated the unstable dif- 
fusive mode. Similarly, the work of Hu & Sugiyama (1996) 
provides a concise analytic treatment of acoustic mode damp- 
ing, but again neglects the importance of gravity at large co- 
moving wavenumber. For the unstable mode I have high- 
lighted, the ^-dependence drops out in the high-A: limit be- 
cause of the ^-dependence of the diffusion timescale. Thus, 
to first order in the limit of large f it is not consistent to ne- 
glect gravity at high-A:. Yamamoto et al. (1998) have discussed 
the same instability, but it is important to emphasize that the 
timescale for the mode's growth is a constant on small scales 
as long as the tight coupling approximation is valid (eq. [22) . 
In this way, the mode's existence is not a manifestation of the 



breakdown of the tight coupling approximation; even at very 
high f , the instability exists with fixed t^a. 

Most importantly, on all scales smaller than the Jeans length 
and at all times before decoupling, the medium is unstable, 
albeit on a long timescale. Thus, it is not formally correct 
to think of the medium on scales below the horizon as stable. 
Although it is dynamically stable, it is unstable on the Kelvin- 
Helmholtz timescale for a radiation-pressure supported self- 
gravitating medium (~ kjc/G). The instability should oper- 
ate at all scales smaller than the Jeans length for which the 
medium is tightly coupled. The physics of the mode is simply 
that self-gravitating radiation pressure supported media radi- 
ate their total internal energy content at the Eddington limit. It 
is qualitatively different than the mode described by Gamow 
(1949) and Field (1971). 

The physics of this mode is already contained in the popu- 
lar codes CMBFAST and CAMB for calculating the growth of 
perturbations in the tight coupling limit. However, precision 
analytic and numerical studies of primordial perturbations in 
density and temperature that neglect gravity on small scales 
may find systematic small differences with respect to more 
complete calculations that are attributable to the growth of 
this mode (e.g., in the transfer function [e.g., Yamamoto et 
al. 1998], or the BAO scale, or in the initial conditions for the 
formation of very small scale dark matter halos). 

As noted in ^2j the instability identified in equation ( 1221 
should operate at essentially all times before decoupling be- 
cause it depends only on the ratio (p m /pb) and the opacity 
kj. However, because the growth timescale is constant, it be- 
comes increasingly long compared to the age of the universe 
at earlier times (see eq. l24l ). Thus, one expects it to be most 
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important at the latest times for which the assumptions apply: 
the largest Jeans-stable modes at a time near decoupling. Even 
here, the growth timescale is much longer than the age of the 
universe, as can be seen from equation (124-b . More generally, 
the medium should be unstable to analogous modes any time 
it is supported against self-gravity by th e rad iation pressure 
provided by a diffusing particle species Q2.21 >. 
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